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2i : Abstract 

Tj • We show that a geometrical notion of entropy, definable in flat space, governs 
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the first quantum correction to the Bekenstein-Hawking black hole entropy. 
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^ ■ We describe two methods for calculating this entropy - a straightforward 

C^* ■ Hamiltonian approach, and a less direct but more powerful Euclidean (heat 

(^ . kernel) method. The entropy diverges in quantum field theory in the absence 

On ■ of an ultraviolet cutoff. Various related finite quantities can be extracted with 



further work. We briefly discuss the corresponding question in string theory. 
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I. CONCEPT OF GEOMETRIC ENTROPY 

The microscopic definition S = — tip In p of entropy, wliere p is tlie density matrix of 
a quantum-mechanical system, has an appeahng information-theoretic meaning apart from 
any thermodynamic interpretation. In quantum field theory, one can define the "geometric 
entropy" associated with a pure state and a geometrical region by forming the pure state 
density matrix, tracing over the field variables inside the region to create an "impure" 
density matrix and then evaluating S. In recent years, 't Hooft [T|] and several other authors 
10-^ have suggested that quantum-mechanical geometric entropy might be related to the 
mysterious Bekenstein-Hawking thermodynamic entropy of a black hole. In this Letter we 
will show that geometric entropy is the first quantum correction to a thermodynamic entropy 
which reduces to the Bekenstein-Hawking entropy in the black hole context. It is therefore 
quite disturbing that, as first noticed by 't Hooft, this quantity has a divergence that cannot 
be renormalized away! 

We begin by presenting a slightly unconventional approach to the Gibbons-Hawking 
black hole entropy. The general expression for thermodynamic entropy is 

S = -(/3^-l)lnZ (1) 

where Z = tre~^^ is the partition function and H is the generator of time translation 
invariance. If the underlying system is a theory of fields {0i(x, t)}, Z can be computed by 
doing the Euclidean functional integral over fields periodic under t ^ t + f3. The periodic 
identification of the time coordinate gives a topology to the Euclidean spacetime: Si x R^, 
for standard thermodynamics and 5*2 x i?2 for the black hole. In the standard case, (3 can be 
given any value and the entropy can be computed for any temperature. Things are not so 
simple in the black hole case: The radial coordinate runs from rhorizon to infinity, while r is 
an angle of period f3. For general (3, there is a two-dimensional conical curvature singularity 
at the horizon of deficit angle 6 = 2'k{i3h — (3)/ I3h and the geometry is non-singular [i.e. 
5 vanishes) only at the Hawking temperature. This is the right physics: a black hole is in 



thermodynamic equilibrium only at a specific temperature. Recasting the general expression 
for entropy (Eq.|l|) in terms of deficit angle rather than /3, one obtains 

Sbh = i27i— + l)\nZ\s=o. (2) 

Since this expression involves derivatives with respect to deficit angle, it evaluates the Eu- 
clidean action on geometries that are not strict solutions of Einstein's equations, and one 
might worry that it is ill-defined. However, Teitelboim f^ and Susskind ^ have shown 
that the classical black hole action can be uniquely extended to conical geometries, thus 
giving precise meaning to Eq. |[ They have further shown that it generates the familiar 
Bekenstein-Hawking entropy when applied to any of the standard black hole solutions. 

Consider now a superficially quite different quantity: the quantum entropy of the ground 
state of a free field, 0(x, t), generated by tracing over fields in the Xi < half space "inside" 
the boundary surface Xi = 0. In the basis where the field itself is diagonal, the ground state 
wave function is given by the path integral 

The fields 0(x, r) are defined on all space and all positive imaginary time, and are subject 
to the boundary condition 0(x, 0) = 0(x). The density matrix p(0^,0^), where 0+ are 
defined for xi > only, is given by a product of two integrals of this type, with an extra 
trace over a common boundary value 0_ on Xi < 0: 

If we let the imaginary time in the two path integrals run from to ±oo, respectively, the 
density matrix becomes the path integral over fields defined on the full Euclidean space, 
with the (xi,r) plane cut along the Xi > axis and with taking values (j)_^ above and 
below the cut. The trace of the density matrix is obtained by equating the fields across the 
cut and doing the unrestricted Euclidean path integral. By the same token, the trace of the 
n-th power of the density matrix is given by the Euclidean functional integral over fields 



defined on an n-sheeted covering of tlie cut spacetime. Note tliat tlie n-fold cover turns tlie 
{xi, t) plane into a flat cone with deficit angle 27r(l — n) at the origin. 

This geometrical interpretation gives a natural definition of the analytic continuation of 
trp" to arbitrary n which in turn allows us to use a "replica trick" to evaluate the entropy. 
In order to eliminate any concerns about normalization, we will construct the entropy out 
of the explicitly normalized object p = p/trp. Then we can write 

'S'geom = -tr(plnp) = (-— + 1) lntrp"|„=i . (5) 

This expression has the useful property of being independent of overall rescalings of p and 
we will, from now on, be rather cavalier about the normalization of p. The identification of 
trp" as the partition function, Zs, on a space of deficit angle 5 = 27r(l — n) finally allows us 
to rewrite Eq.|^ as 

•S'geom = {2n— + l)lnZs\s=o- (6) 

But this is formally the same as the expression Eq.|^ for the black hole thermodynamic 
entropy! To be more precise: the Bekenstein-Hawking entropy and the geometrical entropy 
are the classical and first quantum contributions to a unified object measuring the response 
of the Euclidean path integral to the introduction of a conical singularity in the underlying 
geometry. 

There is one geometric subtlety to mention: In the free field geometric entropy problem, 
the conical geometry obtained by multiply covering fiat space is fiat, except for the coni- 
cal singularity itself. By contrast, the conical geometry obtained by giving the Euclidean 
Schwarzschild solution the wrong f3 has bulk curvature. If the Bekenstein-Hawking and 
geometrical quantum entropies are to be just different orders of approximation to the same 
thing, the conical geometries must be exactly the same. We will achieve this by taking the 
large mass limit, in which the Schwarzschild metric goes into Rindler space. In this limit, 
curvatures go to zero, the area of the horizon goes to infinity (so that the quantity of interest 
becomes entropy per unit horizon area), and the two conical geometries match precisely. Our 



theorem, then, is that the (appropriately defined) geometric entropy of a free field in flat 
space is just the quantum correction to the Bekenstein-Hawking entropy of Rindler space. 
By introducing curvature in the space on which the free field lives, we could compute the 
quantum correction to the entropy of a finite mass black hole. 

As Holzhey and Srednicki discovered, and we shall further discuss below, the geometric 
entropy actually diverges in the absence of an ultraviolet cutoff. Thus the most important 
correlations are those between points just inside and just outside, that is very near to the 
horizon, where flat space is an adequate approximation. The essential physics responsible for 
the diverging quantum correction to the entropy is not some exotic high curvature quantum 
gravity effect, but rather the existence of strong correlations between nearby fleld variables in 
flat space - a basic tenet of generic special relativistic quantum fleld theories. In projecting 
the wave functions onto the ground state as above, we have implicitly assumed regularity 
across the horizon. This means, in the black hole context, that we have chosen the Unruh 
vacuum. 

II. GEOMETRIC ENTROPY BY REPLICAS 

We now turn to a more explicit computation of geometric entropy for free bosonic flelds. 
We will use the replica method directly, in order to give an independent check of the conical 
geometry method described above. By and large, we will be rederiving old (if not widely 
known) results p|, but the exercise will be an instructive illustration of the workings of the 
replica method - and a veriflcation, in our context, of its validity. 

In Eq.0 we wrote the ground state wave functional as a Gaussian functional integral, 
subject to boundary conditions at r = 0. Doing the integral by Green's function methods 
gives the result 



'^o{(f){x)) oc exp{— / (f){x)T{x — x')(j){x')dxdx'} (7) 

where the kernel F is related to the Dirichlet Green's function (satisfying G{x, t = 0;y) = 0) 
by 



1 ^2 

r(a;,y) = -7r^GD(x,r;|/, 0)1^=0 • 



(8) 



2 9r2 

If we split the base space into two regions and correspondingly split the field variables into 
"inside" {(f)-) and "outside" (0+), we can, in an obvious compressed notation, rewrite the 
ground state wave function as 

' 7++ 7+- 

.7-+ 7—/ \0- 
where 7-1--1- denotes the kernel with arguments restricted to "inside" or "outside" in various 

combinations. The integral for the ground state density matrix, Eq.^, can be done explicitly 

giving 



\I'o oc exp{- 
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(9) 
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where 



A = 2(7++-i7+„(7__)-V+) 



B 



h+-{i—) ^7-+ 



(11) 



where the inverse kernels are inverses of the restricted kernels, not restrictions of the full 



Using Eq.|IO|, we can write trp" as a functional integral: 



trp" = f[V^+]exp< 
where A^^ is the supermatrix 
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(13) 



and A and B are the kernels defined in Eq.lTTl Doing the Gaussian functional integral gives 



the result trp" = 1/vdetA^. The problem is to compute this functional determinant as an 
explicit function of n. 

According to the discussion following Eq.^, we can freely rescale p without affecting the 

value of the entropy. Let us define a new matrix 

/ 2 -C ... -C\ 

-C 2 -C ... 

-C 2 ... 
■^n= . . . . . (14) 

... 2 -C 
\-C ... -C 2 / 
with C = —2A^^B. C can be more symmetrically expressed as C = (2 — D)^^D where 



D = 7++7+-7_i7- 



(15) 



is a symmetric matrix acting on the space of 0+ variables. Since detA^^ = 
2'^detA4n/ (detA)^ , replacing J\4n by Ain amounts to rescaling p by a factor of |detA and 
has no effect on the entropy. 

If the operator C were a c- number, we could explicitly diagonalize A^„ and show that 

,27rr. 



det>l„=n{2-2Ccos(— )} = 2"|^^ 



(16) 



where ^(C) is defined by C = 2^/{^'^ + 1). Now C is an operator, not a number, but we may 
diagonalize it and apply the above reasoning eigenvalue by eigenvalue: 



detA^ -2"TT ^^~^^^'^"^^ 



aa 



(17) 



For the usual reasons, we may drop the factors of 2"(1 — C,"^)^ since they have no effect on 
the entropy. The operations called for by the replica trick (Eq.^) then give the entropy as a 
sum over contributions from each eigenvalue of C: 
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(18) 



Although it is not obvious, the eigenvalues of C all lie in (0, 1), as do the values of ^(C), so 
the sum for 5* is positive term-by-term. 

The essential point of all this is that the geometric entropy is calculable from the spectrum 
of a kernel C, acting on the "outside" space. Actually, C is constructed out of the standard 
Green's functions in a very awkward way, and it is better to look at the spectrum of a 
more convenient function of C . The operator D, defined in Eq.|l5|, and related to C by 
C = (2 — D)^^D, is one possibility. It turns out that the operator E = {1 — D)~^D is much 
nicer: it can be expressed as a single convolution 

E ^ -(r-i)+_ r_+ (19) 

of the original kernel F, defined in Eq.^, with its inverse (as usual, the subscripts indicate 
restrictions of the arguments to "inside" or "outside" values). The eigenvalues of E are 
positive and the relation to ^ is 

aE) = i + ^ii-VTTE). (20) 



The entropy sum, Eq.|18|, can of course be expressed directly in terms of the eigenvalues E: 



^replica = Eil^V V + Vl + ^M^ + jl + h} (21) 

Ei ^ V^i V -^« 



III. EXPLICIT CALCULATIONS 

Consider a one- dimensional free field of mass /i, confined to a box —L < x < L to 
eliminate infrared problems. We want to compute the entropy of the ground state density 
matrix obtained by tracing over the fields in —L < x < (so that "inside" ("outside") 
correspond to a; < (x > 0)). We have outlined two methods, Eq.|l^ and Eq.^, which we 
might call Hamiltonian and Euclidean, respectively. They should agree with each other and 
with conformal field theory, where applicable 0. 

Let us first take the Hamiltonian approach and try to compute the eigenvalue sum in 
Eq.|l8|. The most straightforward approach is to discretize the system and then solve the 

8 



resulting finite-dimensional eigenvalue problem numerically. Although this calculation is 
instructive, and we have done it to double-check certain things, it turns out to be hard 
to extract the continuum-limit behavior of 5* from it. Let us instead try an approximate 
analytic approach. In the absence of the infrared cutoff (L -^ oo), the kernels needed to 
construct E are explicit and simple: 

r(a;,?/)"^= Ko{ij,{x-y))/TT 

r(x, y) = fiK,{fi{x - y))/{n{x - y)) (22) 

The eigenvalue problem to be solved is then (this eigenvalue problem appears in [^; our 
contribution is only to show its connection, via the replica trick, with the black hole entropy 
problem) 

Hdy r dzKoiK^ + z)) ^,'^^[' ^/\V (l/) = E^l;ix) (23) 

Jo Jo TT^fi~^[Z + y) 

where a; > 0. It becomes especially simple in the /x — > limit: 

- r dy r dz^^^f^±%^{y) = Eij{x). (24) 

JO Jo T^^[y + zy 

One can easily show by contour integration that ip{x) = exp{tujlnx) is an eigensolution with 
eigenvalue E = sinh~ (vrcj) . The log-periodic nature of these eigenvectors is very reminiscent 
of the Unruh basis for the calculation of the Hawking radiation. 

To calculate the entropy, we must discretize the spectrum. This is normally done by 
imposing boundary conditions. A Dirichlet boundary condition at some large x = L is 
needed, in any case, to regulate infrared divergences. It seems natural, in addition, to 
impose an ultraviolet cutoff via a further Dirichlet condition at some small x = e (this 
is just a reasonable guess: the integral equation picks its own boundary condition!). The 
eigensolutions are then 

iPe{x) = sm{uj{E)\n{x/e)) 

uj{En) ln(L/e) = vrn , (25) 



where u}{E) = vr^^ sinh~ (yE) is the eigenvalue function. It also seems reasonable to keep 
a finite number of degrees of freedom by cutting off the mode number n at L/e. The picture 
so obtained agrees very nicely with the numerical work we have done on the discretized 
problem. The crucial point is that, in the limit L/e —>■ oo, the eigenvalues become dense 
with an eigenvalue density 

-^piE) = -^-^HL/e) (26) 

which diverges logarithmically with the cutoff. 

Note that the density of states per unit u interval is constant, suggesting that it might 
make more sense to rewrite the continuum limit of the entropy expression, Eq.^ as an 
integral over a;, expressing the integrand as a function of u rather than E. The result of 
doing this is 



S„„r..,= -^ d^{-^——-\n{l-e-"")}. (27) 



71 Jo (e^'^'^ - 1) 

an expression which is easily recognized as the entropy of a one-dimensional gas of massless 
bosons confined to a box of size In(-) at a temperature /? = 27r. The integral can be done 
exactly (see any statistical mechanics text) and yields 

^replica ~ ^ln(L/e) + 0(L°). (28) 

The essential fact is that the entropy diverges, not because the contribution of large or 
small eigenvalues blows up, but rather because the eigenvalue density blows up everywhere 
like ln(L/e). This is a direct refiection of the Unruh-style log-periodic behavior of the 
eigenvectors. 

Some further remarks are in order. The formula in Eq.p^ is precisely the expression 



derived by Kabat and Strassler [^] for the black hole entropy in the large-mass (Rindler 
space) limit. This is another piece of evidence for the identity of geometric and black hole 
entropy. Also, Eq.|2^ cries out to be generalized to fermions: the corresponding formula for 
the entropy of massless one-dimensional fermions at zero chemical potential (to pick out the 

10 



vacuum state) is obtained by a strategic sign change reflecting Fermi-Dirac as opposed to 
Bose-Einstein statistics: 



S, 



replica 



In(^) /•«= 27ruj 



r da;{--±!^ + ln(l + e-2--)}. (29) 

n Jo '-(e^'^'^ + l) V ;j \ / 

This integral can also be done exactly, giving the result 

SfZui - ^HL/e) + 0{L'). (30) 

It thus seems likely that the divergence of the geometric entropy will not be cancelled by 
any supersymmetric Bose-Fermi cancellation. A direct calculation of the fermionic entropy 
would be interesting to do. 

Next, let us use heat kernel methods to compute the functional determinant appearing 
in the Euclidean expression, Eq.^ for the entropy. This calculation can easily be done in 
arbitrary dimension: we study a Euclidean fleld theory on the space Cs x Md-2 where Cs is 
a two-dimensional cone of radius L, deflcit angle 6 and Md_2 is a flat {D — 2)-dimensional 
transverse space with total volume Vd-2- For a free massive scalar fleld, the path integral 
is determined by the functional determinant of the kinetic energy operator: 

\nZs = -- lndet(-A + /i^) . (31) 

To calculate the determinant, we use the eigenvalues A„ of —A on the manifold under study 
to deflne the "heat kernel" 

C(t)=tr(e*^)=^e-*^". (32) 

n 

( depends implicitly on the deflcit angle 6 of the geometry. A cut-off version of the deter- 
minant can then be deflned as 

fOO rjf 

In det (-A + f,^) = - -C^{t)e~^' (33) 

J e^ t 

where e is a short distance cutoff with dimensions of length. The heat kernel expression for 
the entropy is then 

11 



The utility of the heat kernel stems from the fact that its expansion in powers of t 
near t = is governed by simple geometrical properties of the manifold |0,^. On the 
two-dimensional manifold Cs, for example, 

a,..-.(!^.^,-L,-_|„0,/iV (35) 

D — 2 

On the full manifold Cg x M£)_2, there is an extra factor VD_2/(47rt)^^ coming from the 
trivial Laplacian on the transverse space. Substituting this into Eq.^ we get an expansion 
of the entropy in powers of L~^: 

= (l/6)ln(l/e) + . . . p = 2) 

= (2V^AVi) ^'"^ + --- (^>2)- (36) 

Note that the replica trick differentiation with respect to deficit angle annihilates the 
^-D/2 ^gj^j^ \^ C,{t). This would have given an extensive entropy proportional to the total 
volume of the D-manifold. The surviving t° term gives a boundary entropy proportional 
to the volume of the boundary manifold My)_i. It is logarithmically divergent in D = 2 
and correspondingly more divergent in higher dimensions. The divergent piece in D = 2 
agrees with our previous "direct" calculation and with |^ while the quadratically divergent 
boundary area entropy in D = 4 agrees approximately with Srednicki's numerical evaluation 
0]. In fact it is unlikely that the coefficient of the leading divergence is universal - that 
is, independent of the regulator procedure - ioY D ^ 2, since it depends on the scale of e. 
Further terms in the expansion of C^ii) give finite volume corrections which vanish in the 
limit L ^ oo. 
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IV. COMMENTS 

1. There are also finite quantities lurking within this circle of ideas which seem likely to 
be of physical significance. One class of such quantities is the difference between geometric 
entropy for the same region but with respect to different states (for example, thermal or 
solitonic states). Another class is the geometric entropy for the ground state in its depen- 
dence on the shape of the bounding surface. Specifically, we suspect that one can formulate 
an expansion for the geometric entropy associated with a bounding surface similar to an ef- 
fective action, where the coefficients of invariants associated with high dimension operators 
{i.e. powers of curvatures and their derivatives) will be cutoff independent. Another arises 
in the dependence on the geometry of the embedding space, for a given bounding surface. 
Yet another occurs in topological field theories, where the geometric entropy associated with 
excision of a topologically non-trivial subspace defines a finite topological invariant. The 
heat kernel algorithm appears to provide a practical approach to such questions. 

2. String theory is thought to provide a consistent framework for calculating quantum 
corrections to classical general relativity, and in particular to be free of the ultraviolet 
divergences that plague standard local field theories in this regard. It is therefore of great 
interest to consider quantum corrections to the Bekenstein-Hawking entropy in the context 
of string theory p. While a complete solution to this problem is far out of sight, our 
considerations here suggest a subproblem that appears to us both central and tractable. As 
we have seen, there is a divergent piece of the entropy arising in order g'^ (free field theory) 
that can be identified even from the partition function in fiat space, or more precisely in 
fiat space with a weak conical singularity. In future work we will attempt to get at this 
by studying orbifolds [l^ formed from fiat space by identifications under discrete rotation 
groups Cn, by analytic continuation in n. 

3. How should one react to the divergence of the first quantum correction to the 
Bekenstein-Hawking entropy? Perhaps one's first instinct is to attempt to subtract or renor- 
malize it away. But subtraction is impossible because of the non-trivial dependence of the 
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correction on area; while renormalization, even if possible, would leave the numerical value 
of the entropy indeterminate. In any case we think that our physical interpretation of the 
corrections as reflecting concrete measurable correlations across the horizon precludes any 
such procedure. Thus we are forced, within quantum field theory, either to abandon the 
equilibrium formula dM = TdS or to let the temperature vanish. In string theory the cor- 
rection to the Hawking result may be finite, but even then it need not be small numerically. 
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